Terahertz Radiation Induced Ballistic Electron Transport in Graphene 
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We investigated ultrafast carrier dynamics in graphene with near-infrared transient absorption 
measurement after intense half-cycle terahertz pulse excitation. The terahertz electric field efficiently 
drives the carriers, inducing large transparency in the near-infrared region. Theoretical calculations 
using the Boltzmann transport equation quantitatively reproduce the experimental findings. This 
good agreement suggests that the intense terahertz field should promote remarkable impact ioniza- 
tion process, which leads to suppression of optical phonon emission and results in efficient carrier 
transport in graphene. 

PACS numbers: 72.80.Vp, 42.65.-k, 72.20.Ht, 72.40.+W 
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Graphene possesses many unique properties arising 
from its truly 2D honeycomb structure, and offers us 
a solid state playground of Dirac fermions[l|. One ex- 
citing property in graphene is its high mobility even at 
room temperature owingto its massless band structure 
and high Fermi velocity [2|, leading to long-range ballis- 
tic transport which is significant for both fundamen- 
tal physics and industrial applications. Although great 
achievements of DC electron transport in graphene @- 
6], current saturation restricts electron transport under 
a high electric field with field strength over 10 kV/cm, 
and has been attributed to opticalphonon emission in 
high-energy state of Dirac electrons [5j-[8[. For anticipated 
graphene based nanoelectronics devices, the phonon scat- 
tering may restrict high field operations. 

Thanks to recent advances in high-power and time- 
resolved terahertz (THz) technologies [9|, we can apply 
electric fields of high strength (MV/cm) with a pulse 
width of ^100 fs to probe the ultrafast nonlinear optical 
properties of materials in the THz frequency region. For 
instance, novel carrier dynamics in semiconductors have 
been reported such as THz-field-driven mass anisotropy 
in doped InGaAs [lfj, and highly efficient carrier multi- 
plication in GaAs quantum wells [11], to name only two. 
In graphene, S. Winner 1 et al. has performed THz-pump 
THz-probe measurement, revealing that the excitation 
and relaxation dynamics of carriers strongly depend on 
the excitation frequency in the THz region [12]. In these 
studies, THz pulses are of picosecond duration, which is 
longer than the typical optical phonon emission timescale 
in graphene. Shorter pulse is required to prevail over cur- 
rent saturation caused by optical phonon emission. 

In this work, we succeeded in ballistic electron trans- 
port in graphene using 200-fs half-cycle THz excitation 
pulses and monitoring their dynamics with 50-fs long 
near-infrared (NIR) probing pulses. Large THz induced 
transparency is observed in the NIR region for the first 
time. Numerical calculations using the Boltzmann equa- 



tion quantitatively reproduce the experimental findings 
and demonstrate the importance of efficient carrier mul- 
tiplication by impact ionization in graphene. This mul- 
tiplication achieved by ultrashort THz pulse excitation 
drastically suppresses the current saturation effect. 

The sample used in our experiment is a commercial 
CVD graphene grown on Cu foil, which is transferred 
onto a Si02 substrate. We have confirmed the sample's 
single-layer nature and uniformity by confocal Raman 
imaging microscope (l3j . The transmissivity of our sam- 
ple at 800 nm is 97 %, comparable with the value ob- 
tained from universal conductivity of intrinsic single layer 
graphene [l[. Because graphene is deposited on Si02, 
substrate-induced inhomogeneity at the graphene-oxide 
interface gives jdse to p-type doping with Fermi energy 
over 200 meV 
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To investigate hot electron transport properties in 
graphene, we performed transient absorption measure- 
ment in the NIR region under intense THz excitation, as 
illustrated in Fig. HJa). We used a collinear air plasma 
method for generation of intense half-cycle THz pulses 
[3 (Fig. QJb)). A 800-nm laser pulse (Ti: Sapphire laser 
that delivers 3 mJ with 35-fs duration laser pulses) was 
split into two pulses for THz excitation and time-resolved 
NIR transmissivity measurements. An undoped germa- 
nium substrate was inserted in the THz excitation line 
to reject any remaining 800 nm-pump beam after plasma 
generation. The laser pump and probe energy are 2.4 mJ 
and <1 nJ per pulse, and spot sizes of THz pulse and 800- 
nm pulse at the sample position are 80 /^m and 30 /im, 
respectively. Figure QJc) shows a temporal profile of the 
THz pulse with a peak electric field of 300 kV/cm and a 
full-width at half-maximum of ~ 100 fs. Notice that spec- 
trum of THz pulse is ranging from 0.5 to 7 THz (inset of 
Fig. [He)) w hh the photon energy several tens of times 
smaller than that of the probing pulse. All measure- 
ments were performed at room temperature and under 
dehydrated air to suppress THz absorption from water 
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FIG. 1. (a) Schematic figure of THz-pump NIR-transient ab- 
sorption measurement in Dirac cone, (b) Experimental setup, 
(c) Incident THz electric field and its spectrum, (d) THz in- 
duced normalized differential optical density of graphene at 
800 nm as a function of delay time. 



vapor. Temporal evolutions of 800-nm probe transmis- 
sion through graphene after THz pulse excitation could 
be reconstructed by scanning an optical delay stage. Fi- 
nally, a lock-in amplifier (SRS SR830) connected to the 
output of an amplified photodiode, and referenced to the 
chopper that switched the excitation on and off at 500 
Hz, served to probe the change of the NIR transmissivity 
AT. 

Figure QJd) shows a typical THz induced differ- 
ential optical density measurement as a function of 
time with a peak THz electric field of 300 kV/cm. 
To clarify the physical meaning observed here, we 
define the normalized differential optical density by 
AOD/OD rof = (OD ox - OD re f)/OD re f ~ AT/OD rof , 
where OD ro f is the optical density of graphene at 800 
nm without THz pulse excitation, and OD cx is that of 
with THz pulse excitation. The value of — AOD/OD re f 
corresponds to the population occupation at the energy 
of 800-nm transition. After THz pulse excitation, the dif- 
ferential optical density immediately goes negative (e.g., 
graphene becomes transparent in the NfR region un- 
der THz pulse excitation) and returns to zero within 2 
ps. Experimentally, the maximum value of AT reaches 
5 x 10~ 3 for a peak electric field of 300 kV/cm. This large 
induced transparency over 16% was probably caused by 
to electron filling or hole depletion of the corresponding 
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FIG. 2. THz induced normalized differential optical density 
as a function of time delay with peak electric fields of (a) 170 
kV/cm and (b) 300 kV/cm : experiment (circle) and sim- 
ulation (solid line). (c)(d) Calculated carrier density as a 
function of time delay. Dashed lines represent initial carrier 
density, (e) THz induced transparency of graphene as a func- 
tion of a peak electric field : experiment (filled circle) and 
simulations (solid and dashed lines). Experimental results 
follow a power law with exponent 2.5 (dotted line). 



energy levels. 

Figures (Ha) and HJb) show the temporal profile of in- 
duced NIR transparency for peak electric fields of 170 
kV/cm and 300 kV/cm. The maximum values of induced 
transparency drastically increase for the higher electric 
field, while carrier relaxation time tends to increase lin- 
early with field strength. The origins leading to the non- 
linearity will be discussed in subsequent paragraphs. 

Figure [5Je) shows the maximum values of induced 
transparency as a function of a peak electric field. A 
power law with exponent 2.5 is found experimentally. 
This is in contrast with quadratic dependence observed in 



ultrafast optical pump-probe experiments |16l |17| . The 
super-quadratic dependence found here suggests nonlin- 
ear coupling between the THz field and electronic system, 
and can be explained by the rapid increase in the number 
of carriers with field strength. 

As examined in detail in previous studies from THz to 
NIR regions 12 , 1|| 1{| , electromagnetic wave excitation 
results in two possible scenarios: interband and intra- 
band transitions. Interband transition creates electron- 
hole pairs with kinetic energy of the incident photon en- 
ergy, whereas intraband transition induces athermal dis- 
tribution of carriers in the momentum space. Accord- 
ing to these studies, the dominant mechanism should 
be determined by photon energy and Fermi energy of 
graphene. In this study, the dominant mechanism is re- 
stricted to intraband transition because of large nega- 
tive Fermi energy below -200 meV coming from the SiC>2 
substrate. Since the excited carriers have the kinetic en- 
ergy corresponding to NIR excitation, their relaxation 
dynamics can be compared with photo-excitation studies 
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FIG. 3. (a)-(e) Distribution function of holes in k-space with various time delays. Dashed lines represent momenta corresponding 
to the energy of 800-nm transition, (f)-(j) Angle-resolved energy distribution of holes. Fermi-Dirac fitting is represented as a 
dashed line. Brown arrows indicate the energy of 800-nm transition. 



such as ultrafast saturation absorption [20f. sin gle-c olor 
and two-color pump-probe measurement [lot Il7|. l2l H23 |. 
and ultrafast photoluminescence measurement [24| , where 
8 fs carrier-carrier scattering, 50-120 fs carrier thermal- 
ization, 0.1-1 ps carrier cooling and 3 - 15 ps carrier re- 
combination have been reported . 

To understand the carrier dynamics based on our ob- 
servations, we have calculated the distribution function 
of carriers in k-space using the following semi-classical 
Boltzmann equation: 



|A(M) 



the second term, momentum-resolved hot phonon effects 
are included 26[. Acoustic phonon scattering negligibly 
contributes in our observation time scale and was not 
included in our calculations (l2. 26 1. 

The occurrence of II and AR in graphene is controver- 
sial because of restrictions from energy and momentum 
conservation laws in the Dirac cone. Although some au- 
thors have stated that these processes are prohibited be- 
cause of a limited number of possible initial and final mo- 
mentum configurations [2^1, some experimental and the- 
oretical studies have indicated the possibility of using 
these processes [IT 
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t these processes by loosening the momentum conserva- 

k',k ~ ^k',k } tion law by Afc = 2 x 10 4 cm -1 (~ 1 meV in kinetic 
energy) . This may come from imperfection of our CVD 
graphene, which possibly results in larger II and AR pro- 
fit " 



cesses 
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where f\(k,t) is a distribution function of electrons and 
holes (A = e,h), q\ is the charge of carrier, and E{t) is the 
incident THz electric field on graphene at time t. S k ~^, 
represents intraband carrier-carrier (c-c) scattering rates 
with coulomb interactions [25j. Intervalley K-K' scat- 
tering is neglected for large momentum exchange [2fjj |. 
S k ^? represents intra- and inter-band phonon scattering 
rate with T and K optical phonon modes [26]. Electron- 
phonon coupling strength is from Ref. |27| . SfJ£ and 

Sufu^ represent impact ionization (II) and Auger recom- 
bination ( AR) rates |28[ , respectively. The last two terms 
(II and AR) change the number of carriers through in- 
terband c-c scattering. The II process creates an extra 
electron-hole pair while losing the kinetic energy of an- 
other electron or hole, as schematically shown in Fig. 
[3j The AR is the inverse process of II, which reduces 
the number of carriers. The first term represents car- 
rier acceleration on linear dispersion relation E — hvp\k\ 
(Dirac cone) with an electric field, where E is the energy 
of the carrier and vf is the Fermi energy in graphene. In 



Calculated induced NIR transparency is plotted as a 
function of a peak electric field in Fig. [2je) (solid line). 
The experimental results are quantitatively reproduced 
when II and AR processes are included. For comparison, 
the calculations by strictly applying conservation laws 
(without II and AR) are also plotted as a dashed line. In 
the latter case, the sub-linear scaling with field strength 
clearly deviates from the experiment for the higher elec- 
tric field because of rapid optical phonon scattering. We 
have to stress here that the only adjustable parameter 
is the initial Fermi energy, which is tuned to reproduce 
the magnitude of induced transparency in Fig. (He). The 
best-fitted Fermi energy is -280 meV, which is a reason- 
able value for graphene on a Si02 substrate |l4[. Note 
that the changes of the initial Fermi energy only offset the 
field dependence without affecting the power law. This 
good agreement indicated that efficient carrier multipli- 
cation by II prevails over AR under the high electric field 
condition. 

To emphasize the effect of carrier multiplication, we 
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plot the calculated carrier density as a function of time 
in Figs. [2c) and[2jd). The number of carriers becomes 
almost five times larger than the initial number of car- 
riers with a peak electric field of 300 kV/cm, resulting 
in the observed large induced transparency. Recovery of 
the absorption change at 800 nm is also calculated and 
shown as solid lines in Figs. [2ta) and[UJb). Our calcu- 
lations qualitatively reproduce the field dependent relax- 
ation time, which confirm the longer relaxation time at 
higher excitation density. This is due to the hot phonon 
effect 0,[26|. 



Figures E3a)-(e) show calculated distribution of holes 
in k-space at various time delays for the maximum ap- 
plied electric field of 300 kV/cm. Figures Et^f)-(j) show 
the corresponding angle-resolved energy distribution of 
holes defined as fg(E) = 0))0_ 7r /4<<£<0+7r/4 with 

the angle 9 in Fig. Efa). Before THz pulse excitation, 
holes are distributed within Fermi energy with 300 K 
thermal fluctuation (Figs. |3fa) and EIO)- As a THz 
electric field is applied, the carrier distribution becomes 
an asymmetric shape (Figs. [3Jb) and^g)), and keeps its 
asymmetry until the peak electric field (Figs. (3Jc) and 
[31(h)), although the asymmetry decreases rapidly by op- 
tical phonon scattering. After the THz excitation pulse 
is passing through the sample, the distribution immedi- 
ately recovers its symmetric shape (Figs. EJd) and[3ji)) 
and returns back to the initial distribution by emitting 
optical phonons (Figs. Efe) and[3Jj)). This argument 
is well confirmed by angle-resolved energy distributions 
fitted with Fermi-Dirac distribution in the THz field di- 
rection k x (Figs. [3jf)-(j)). For instance, temperature in 
k x direction is 300 K before excitation and is increased 
to 2860 K after THz pulse excitation, though temper- 
ature in other direction remains low even at the peak 
electric field (Fig. GHh)). After THz pulse excitation, 
energy distributions immediately thermalize in all direc- 
tions and return back to the initial distribution through 
phonon emission. Notice that the dashed circles in Figs. 
ini^a)-(e) represent momenta corresponding to the energy 
of 800-nm transition, also indicated with arrows in Figs. 

The characteristic feature of half-cycle THz pulse exci- 
tation is field-induced asymmetry in the momentum dis- 
tribution, which induces current in the direction of the 
THz electric field. In the linear dispersion relation, a car- 
rier with the momentum direction 9 induces current by 
q\VF cos(#). With this expression, we can estimate how 
the scattering processes affect current. In short, optical 
phonon scattering efficiently damps current, whereas c-c 
scattering, II and AR processes conserve current. 

In contrast to high DC electric field studies in which 
rapid optical emission restricts current flow [H-IH > current 
is not saturated under intense THz pulse excitation. This 
can be understood as follows: In our high field conditions, 
carriers can gain enough kinetic energy to induce im- 



pact ionization before optical phonon emission. Because 
impact ionization process efficiently reduces kinetic en- 
ergies of carriers while conserving current, carriers have 
less chance to emit optical phonons. This suppresses the 
current saturation and leads to longer ballistic transport 
of carriers in graphene under high electric fields. 

In conclusion, we reported the first observation of large 
THz-induced near-infrared transparency in graphene. 
Our results were compared with simulations using the 
Boltzmann equation and we proved that carrier multi- 
plication occurs in graphene driven by high field and ul- 
trashort THz pulse. The observed dependence of the 
rise in charge density on THz peak field strength in- 
dicates an interband impact ionization mechanism in 
graphene. The impact ionization process may suppress 
optical phonon emission and lead to ballistic carrier 
transport in graphene. 
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